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ON THE COMPUTATION OF FACTORIZATION INVARIANTS FOR 

AFFINE SEMIGROUPS 

PEDRO GARCIA-SANCHEZ, CHRISTOPHER O’NEILL, AND GAUTAM WEBB 


Abstract. We present several new algorithms for computing factorization invariant values 
over affine semigroups. In particular, we give (i) the first known algorithm to compute 
the delta set of any affine semigroup, (ii) an improved method of computing the tame 
degree of an affine semigroup, and (iii) a dynamic algorithm to compute catenary degrees 
of affine semigroup elements. Our algorithms rely on theoretical results from combinatorial 
commutative algebra involving Grobner bases, Hilbert bases, and other standard techniques. 
Implementation in the computer algebra system GAP is discussed. 


1. Introduction 


Let r = (cti,... ,afc) C be a finitely generated subsemigroup of that is, an affine 
semigroup. We are interested in studying the factorization structure of elements 7 G L, 
that is, the distinct expressions of 7 as a sum of the generators ai,... ,ak. Most of the 
theory on nonunique factorization stems from the study of factorizations in integral domains; 
in this context, one can view factorizations of affine semigroup elements as factorizations 
of monomials in the corresponding monoid algebra IK[r] for K a field. Often, factorization 
structure is examined using factorization invariants, which assign quantities to each element 
of r (or to r as a whole) that measure the failure of its factorizations to be unique [22] . 

Every affine semigroup T is an FF-monoid, meaning that each element 7 G L has only 
finitely many factorizations. As such, many factorization invariants can be computed explic¬ 
itly using computer software |20|, an attribute that has greatly aided their study |28|. On the 
other hand, the class of affine semigroups is sufficiently broad to exhibit varied factorization 
structure, and remains an interesting source of examples and an active area of study. 

In this paper, we present new algorithms for computing several well-studied factorization 
invariants. After some preliminary definitions in Section]^ we present in Section]^ the first 
known algorithm for computing the delta set of an affine semigroup. Two distinct approaches 
are given (Algorithms |3 . 6 | and |3 . 9[ ) , the latter of which introduces a novel connection between 
the delta set and Grobner bases of homogeneous toric ideals (see Remark 3.12). Next, we give 
Algorithm |4.11| for dynamically computing catenary degrees of affine semigroup elements, in 
the style of several known algorithms for other factorization invariants [3] . We conclude with 
an improved algorithm for computing the tame degree of an affine semigroup in Section 

Following the statement and proof of correctness of each algorithm, we discuss implemen¬ 
tation in the GAP m package numericalsgps [Hj, including sample code and benchmark 
comparisons with existing algorithms. Whenever a Grobner basis (Definition |2.9| ) or a Hilbert 
basis (Definition 2.3[ ) must be computed, our implementations rely on trusted software pack¬ 
ages like 4ti2 [T], Normaliz |5|, and Singular |13j . 
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Notation. In what follows, let N = {0,1,2,...} denote the set of non-negative integers. Fix 
a field K, and let ]K[yi,..., yk] denote the polynomial ring in commuting variables yi, ■ ■ ■ ,yk 
with coefficients inK. Lastly, for z G we use the shorthand y^ = yl^ ■ ■ ■ y^'^ G ^[ui, ■ ■ ■ ,yk]- 

2. Background 

Definition 2.1. A semigroup F C N'^ is affine if it is finitely generated. If F C N and 
gcd(r) = 1, we say F is a numerical semigroup. 

Remark 2.2. Any affine semigroup F is reduced, and thus has a unique generating set 
that is minimal with respect to containment. Throughout this paper, whenever we write 
F = (oi,..., Oik): we assume oi,..., are precisely the minimal generators of F. 

Definition 2.3. Fix a matrix A G The set 

F = {x G : Ax = 0} 

of nonnegative integer homogeneous solutions of A forms an affine semigroup. The minimal 
generating set of F (w.r.t. containment) is called the Hilbert basis of A, denoted 

Definition 2.4. Fix an affine semigroup F = {ai ,..., ak) C N'^. The elements oi,..., 
comprising the unique minimal generating set of F are called irreducible elements (or atoms). 
A factorization of 7 G F is an expression 

7 = zioi H-h Zkak 

of 7 as a finite sum of atoms, which we denote by the fc-tuple z = (zi,..., Zk) G N^. Write 
Zr( 7 ) for the set of factorizations of 7 G F, viewed as a subset of N^. 

Since ai,... ,ak generate F, the monoid homomorphism (/?r ^ F given by 

^prizi, ...,Zk) = zioi H-h ZkOk 

is surjective, and F = N^/kert/ppj where 

ker(/?r = {iz,w) G x : cpriz) = (/ 2 r(R')} 

denotes the kernel of (pr. Notice that ker pr is a congruence on (i.e. an equivalence relation 
that is closed under translation), and thus it is finitely generated (see [29]). In particular, 
there exists a finite set p C ker <^r such that the smallest congruence containing p is ker ip-p 
(equivalently, ker equals the intersection of all congruences containing p). 

Definition 2.5. A minimal presentation of F is a generating set p of the congruence ker (/?r 
that is minimal with respect to containment. The Betti elements of F are elements of the 
form Pt{z) for {z,w) G p (this is independent of the minimal presentation chosen). 

Definition 2.6. Fix F = (oi,..., o^) C N'^, and let A = (oi | • • • | Ofc) denote the matrix 
with columns oi,... ,ak. The Graver basis of A (or the Graver basis of F) is the Hilbert 
basis of the matrix {A \ —A). 
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Remark 2.7. The Graver basis of a semigroup T contains every minimal presentation of T, 
but is typically orders of magnitude larger. Even for numerical semigroups with 3 minimal 
generators, whose minimal presentations have at most 3 relations, the Graver basis can be 
arbitrarily large. As such, algorithms relying on a proper subset of the Graver basis (such as 
a minimal presentation or a Grobner basis) are preferred whenever possible. 

Example 2.8. Let T = (3,4,5) C N. The minimal presentation of T is 

p = {((1, 0,1), (0, 2, 0)), ((2,1,0), (0,0, 2)), ((3, 0, 0), (0,1,1))}, 

as computed in m Example 8.23]. In this case, p is the only minimal presentation of T, 
though minimal presentations need not be unique in general. For comparison, the Graver 
basis of T consists of (upon omitting symmetry) the pairs 

((1, 0,1), (0, 2, 0)), ((1,3,0), (0,0, 3)), ((2,1, 0), (0,0, 2)), 

((3, 0, 0), (0,1,1)), ((4,0,0), (0, 3,0)), ((5, 0, 0), (0,0, 3)). 

We associate to T an ideal Ir C ^[yi,... ,yk] defined as follows. Consider the polynomial 
ring IC[ti ,... ,td\, and more specifically the subring ]K[r] = 0 ..j,grIii this setting, pr has 
a polynomial ring counterpart, namely the ring homomorphism 

ipr : IK[yi,... ,yfc] -A K[T] 

determined by '4’riyi) = t°‘\ The defining ideal of T is the kernel 

Ir = kerfir C K[yi,.. .,yk]. 

It is well known that a set p C ker pr generates ker pr as a congruence if and only if 

{y^ - y^ : (z, re) E p} C h 

generates Ir, and that p is a minimal presentation for T precisely when the above set irre- 
dundantly generates Ir (see [24]L 

We conclude this section with a brief discussion of Grobner bases, which are standard in 
computational commutative algebra; the unfamiliar reader is encouraged to consult [12] for 
a more thorough overview. 

Definition 2.9. Fix a monomial order ^ on IK[pi,... ,pfc] (that is, a total ordering of the 
monomials containing the usual partial ordering by divisibility). Given / E ]K[?/i, ..., p^j, 
denote by In^(/) the initial term of / with respect to A Grobner basis of an ideal 
I C ]K[yi,..., pfc] with respect to ^ is a set G = {gi,...,gr} generating I such that 
(In^(pi),... ,In^(pr)) = ({In^(/): / E /}). 

Remark 2.10. Minimal presentations can be computed using Grobner bases. Set 

T = (pi - C K[yi, ...,yk,ti,... ,td], 

and observe that Ir = Jn]K[pi,..., p^j. To compute a generating set for Ir, one can compute 
a Grobner basis G for J with respect to any monomial ordering ^ satisfying p* ^ tj for all 
possible i and j, and then calculate G n IK[pi,..., p^j. 
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3. The delta set 


In this section, we introduce Algorithms 3.6 and |3.9[ the first known algorithms for com¬ 


puting the delta set of any affine semigroup, as well as proofs of correctness (Theorems 3.8 


and 3.10 respectively). Both algorithms have Theorem 3.5 at their core, but each uses a 


different method for computing the necessary ideal generating sets. 

Definition 3.1. Fix an affine semigroup F = (oi,..., au) C and 7 G F. Given 2 ; G Zr( 7 ), 
the length of 2 ; is the number \z\ = zi + ■ ■ ■ + Zk of irreducibles in 2 :. The length set of 7 is 

■-r(7) = {zi + Z 2 -\ -h 2 ;fc : 2 ; G Zr( 7 )}, 

the set of factorization lengths. Writing Lr( 7 ) = {(-i < ■ ■ ■ < Im}-, the delta set of 7 is the set 

A( 7 ) = {^i+l - : 1 < f < m} 

of successive differences of factorization lengths. The delta set o/F is A(F) = U 7 er^('T)- 


Theorem 3.2 follows from [HI Corollary 2.4, Thoerem 2.5], and describes which delta set 
elements can be easily recovered from a minimal presentation. 

Theorem 3.2. //F = (oi,. .., C and p C is a minimal presentation o/F, then 

minA(F) = gcd({| 2 ;| — |r(;| : {z,w) G p}), 
maxA(F) = maxjmax A( 2 ;iq;i + • • • + 2 :^ 0 ^) | {z,w) G p} 

and 

A(F) C [minA(F),max A(F)] n minA(F)Z. 

In particular, min A(F) = gcd A(F), and both min A(F) and maxA(F) can be quickly recov¬ 
ered from p. 


Example 3.3. Let us return to F = (3,4,5) from Example 2.8 According to Theorem |3.2[ 
the minimum and maximum of A(F) are determined by A(8), A(9) and A(10). The only 
factorizations of 8 are (0, 2, 0) and (1,0,1), so A(8) = 0. We can obtain A(9) = A(10) = {1} 
similarly. As such, A(F) = {1}, since in this case min A(F) = max A(F). In general, however. 
Theorem 3.2 only yields an interval containing A(F). 


Remark 3.4. The primary difficulty in computing A(F) for a general affine semigroup F 
is ensuring that a given value does not occur in A(F). Indeed, some elements of A(F) may 
only occur in the delta sets of a small (finite) number of semigroup elements. For example, 
if F = (17,33,53,71) C N, then 6 G A(F) only occurs in A(266), A(283), and A(300). 

As such, although it is computationally feasible to compute the delta set of any single 
element of F (since each has only finitely many factorizations), this cannot be accomplished 
for all of the (inhnitely many) elements of F. To date, all existing delta set algorithms are 
limited to numerical semigroups, and act by bounding the semigroup elements for which 
“new” elements of A(F) can occur, thus restricting computation to a finite list of semigroup 
elements; see |3] for more detail. 


One of the primary benehts of Algorithm 3.6 is that it does not rely on computing delta 


sets of individual semigroup elements. As such, both of our algorithms are generally much 
faster than existing algorithms for numerical semigroups; see Examples [3.13 and 3.14 














ON THE COMPUTATION OF FACTORIZATION INVARIANTS FOR AFFINE SEMIGROUPS 


5 


We now state Theorem |3.5| the main theoretical result used in Algorithm |3.61 
Theorem 3.5 ( [271 Theorem 4.8]). Suppose T = (ai,... ,ak) C N'^. The ideals 

Ij = {y'' - : z,w G Zr(7),7 e T, and ||u;| - | 2 ;|| < j) C K[yi, ...,yk] 

for j >0 form an ascending chain 


lochchc--- 


in which Ij-i C Ij if and only if j G A(r). 


Algorithm 3.6 obtains the delta set of a given affine semigroup T = {ai,... ,ak) C 


by computing each ideal defined in Theorem 13.51 The generating set of T appearing in 


Theorem 3.5 is determined by (integer solutions to) a system of linear inequalities, and thus 


can be computed using a Hilbert basis. For example, by inserting a slack variable to transform 
the single inequality to an equality, it suffices to compute a Hilbert basis for the matrix A 
below. Each element x = {z, w,i) G x x N of the Hilbert basis with i < j yields a 


binomial generator y^ — y"^ for L by Lemma 3.7 


( ail ■ 

• aifc 

— Oil • • 

• —dik 

0 \ 


• OLdk 

— • • 

■ —Oidk 

0 

V 1 • 

■ 1 

-1 

■ -1 

-1 ) 


A = 


Once each generating set has been computed, successive ideals Ij-i C Ij can be checked 
for strict containment by computing a reduced Grobner basis for each ideal (with respect 
to the same term order) and then comparing. The final step in the algorithm is locating 
the maximal value of j, which by Theorem |3.2| can be obtained by computing a minimal 
presentation for T. Note that each element of A(r) is a multiple of minA(r), which by 


Theorem 3.2 can also be obtained from a minimal presentation for T. 


We are now ready to state Algorithm 3.6 


Algorithm 3.6. Computes the delta set of an affine semigroup T = (oi,..., a^) C 
function DELTASETOFAFFINESEMIGROUP(r) 
p G- minimal presentation for F 

^ ^ U(., ^■^^p/\{ziai H-h Zkak) 

m G- gcd D 

H G- Hilbert basis for {x G : Ax = 0} 

G G- Grobner basis for — y"^ : (z, w,m) G II or (z, w, 0) G H) C IK.[yi ,... ,yk] 
for all j = 2 m,..., maxT) — m do 

if y^ — y^ has nonzero remainder modulo G for some (z,w,j} G H then 
D ^ DVI {j} 

G G- Grobner basis for (G) + (y^ — y"' : (z, w^j) G H) 

end if 
end for 
return D 
end function 


We now give a proof that Algorithm |3.6| gives the correct output. 
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Lemma 3.7. Resume notation from Theorem 3.5 and Algorithm \3.(^ The ideal 
Jj = : {z,w,i) G H,i < j) C ]K[yi, ...,yk] 


computed in Algorithm 3.6 coincides with the ideal L from Theorem\3.^for all j > 0 


Proof. For each i < j and {z, w, i) G H., it is clear that y^ — y^ G Ij, so Jj C Ij. Conversely, 
if (z, w, i), {z', w', i') G {H) with y^ — y^, y^' — y'^' G Ij, then 

y^'{y^ - yn + y^iy^' - y"') = e ij, 

which is the binomial corresponding to {z,w,i) + {z',w',i') G (H) C means the 

generating set for Ij given in Theorem 3.5 can be restricted to the binomials corresponding 
to elements of the Hilbert basis H, so Ij C Jj as well. □ 


Theorem 3.8. For any affine semigroup T C N“, Algorithm 3.6 returns A(r). 


3.2| implies that 


3.7 ensures each 


Proof. Resume notation from Theorem 3.5 and Algorithm |3.6[ Theorem 
m = minA(r), and that only multiples of m can appear in A(r). Lemma 
loop iteration computes the ideal Ij for some multiple j of m, and appends j to D if Ij-i £ Ij- 
Applying Theorem 3.5 completes the proof. □ 


Algorithm 3.9 uses an alternative approach to compute generators for the ideals Ij in 
Theorem 3.5 Given a G denote by (l,a) the element in whose first coordinate is 1 
and whose remaining coordinates are those in a. For F = (ai,..., a^) C N'^, set 


= ((1,0), (l,ai),. 

There is a tight relation between factorizations in T// and those in F 


Observe that 


Irj^ = -y^ -. z,w G ^ T, kl < kl) C K[t,yi,.. .,yk]. 


Algorithm 3.9 exploits this idea to compute the generators L in Theorem 13.51 As there is 


no need to compute a Graver basis (see Remark 2.7), Algorithm 3.9 is better equipped than 
Algorithm 3.6 to handle input with large generators (see Example |3T3l ). 


Algorithm 3.9. Computes the delta set of an affine semigroup F = (oi,..., a^) C 
function DeltaSetOfAffineSemigroup(F) 

p G- minimal presentation for F/^ = ((1, 0), (1, oi),..., (1, Ok)) 

G G- reduced lex Grobner basis for {Py^ — Py'^ : ((i, z), (j, w)) G p) C K[t, yi,... ,yk] 
return {j : Py^ — y'^ G G} 

end function 


Theorem 3.10. Fix F = {ai,... ,ak) C an affine semigroup. Let G denote a reduced 
Grobner basis for with respect to any lexicographic term order satisfying t > yi for eaeh 
i < k, and for j > 0 let 

Jj = {y^-y-’-.tY-y^eG,i<j). 

We have Ij = Jj for each j > 0. 
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Proof. Each generator of Jj certainly lies in Ij , so Jj C Ij . Conversely, fix € Ij with 

\z\ < lu)!, and let i = \w\ — \z\ < j. By definition, Py^ — y"' G Ir^, so performing polynomial 
long division by G yields a remainder of 0. Since G was computed using a lex term order, 
the only elements of G used in the division algorithm have leading term dividing Py^, and 
thus correspond to generators of Jj. As such, Ij C Jj. The final claim follows from G being 
reduced, since the initial term Py^ of each Py^ — y'^ ^ G cannot be divisible by the leading 
term of any other elements of G, including those corresponding to generators of Jj-i. □ 

Corollary 3.11. Fix T = {ai, ..., C N'^ an affine semigroup. If G is a reduced Grobner 
basis for with respect to any lexicographic term order satisfying t>yi for i < k, then 


A{r) = {j-.Py^-y^GG}. 

Remark 3.12. In addition to ensuring the correctness of Algorithm |3.9[ Corollary 3.11 is also 
interesting from a theoretical perspective, as it expresses the delta set in terms of Grobner 
bases of homogeneous toric ideals. The authors are optimistic that this novel connection will 
spur interest in the delta set from members of the computational algebra community. 


3.1. Implementation and benchmarks. We implemented both algorithms above in GAP 
to compare their performance to existing methods. In our implementation of Algorithm |3.6t 
we used Normaliz [5] through Normalizinterf ace |23) for the calculation of the Hilbert basis 
H, and Singular [13] via the GAP package singular [11] for the computation of successive 
Grobner bases. Although we did some experiments with 4ti2 [T] via [2] for the calculation 
of H, the resulting implementation was slower than with normaliz; additionally, one cannot 
use 4ti2 for Grobner basis computations here since the ideals are not toric. 


Example 3.13. For numerical semigroups, the algorithm used in the numericalsgps m 
package is the one presented in [3], which was the fastest known prior to those given above. 
Though it only works with numerical semigroups, we can still compare how it behaves against 
Algorithm |3.6[ A log-time algorithm has also been given for numerical semigroups with 
3 minimal generators |18j . but the techniques used therein do not easily generalize. 

In examples with “small” generators, Algorithm 3.6 outperforms DeltaSetOfNumerical- 
Semigroup (the current implementation in numericalsgps). For instance, the calculation 
of A((101,123,147,199)) takes 550 ms with the current implementation, compared to 50 ms 
with Algorithm |3.6[ Additionally, the current implementation takes 415263 ms to compute 


A((101,123,147,199, 210)), while Algorithm 3.6 accomplishes this task in under 4 seconds. 

For semigroups with large generators, the situation can differ. Ghoosing, for instance, 
T = (1001,1211,1421,1631,2841), DeltaSetOfNumericalSemigroup computes A(r) in less 
than two minutes, whereas the processes running Normaliz completed Algorithm [3.6| in 
around three minutes, even when run with SCIP |15j (4ti2 was killed after half an hour). 

One of the advantages of Algorithm [3^ is that one can use either Singular (eliminate ()) 
or 4ti2 (groebnerO) to compute a Grobner basis for Ipjij, with comparable performance. 

were computed with Algorithm|3.9[ 


Example 3.14. All delta sets appearing in Example 3.13 


in less than 20 ms. Even for significantly larger examples such as 


T = (550,1060,1600,1781,4126,4139,4407,5167, 6073, 6079, 6169, 7097, 7602,8782,8872), 
A(r) = {1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,17,19} computes in less than a minute. 
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4. The catenary degree: a dynamic algorithm 


The catenary degree is a factorization invariant that measures how spread out the factor¬ 
izations of an element are. Recent results prove the catenary degree is eventually periodic for 
numerical semigroups m, and a natural generalization is given for affine semigroups EZj, 
but neither result gives a concrete bound for the start of periodic behavior (in contrast to 
similar results for several other invariants; see |28[ Remark 5.11]). In this section, we give Al¬ 
gorithm |4.11 for computing the catenary degree of affine semigroup elements using dynamic 
programming, answering [281 Problem 5.12]. The primary motivation for such an algorithm is 
to aid the investigation of concrete bounds for the aforementioned eventual behavior results, 
and our approach mirrors several existing algorithms for other factorization invariants [3]. 

Definition 4.1. Fix F = (ai,...,afc) C and 7 G F. The greatest common divisor of 
z,w £ Z( 7 ) is given by 

gcd{z,w) = {mm{zi,wi),... ,m.m{zk,Wk)) G 
and the distance between 2 ; and w (or the weight of {z,w)) is defined as 
d{z, w) = maxdz — gcd( 2 ;, u;)|, |rc — gcd( 2 ;, w)\). 

Given z,w £ Z(j) and > 0, an N-chain from 2 ; to tc is a sequence 2 : 0 ,... , 2 ;^ G Z( 7 ) of 
factorizations of 7 such that (i) zq = z, (ii) Zr = w, and (hi) d(zi-i,Zi) < N for all i < r. 
The weight of 2 : 0 ,..., 2 ;r is the smallest N such that 2 : 0 ,..., 2 ;^ is an A^-chain (or, equivalently, 
the largest distance between successive factorizations). The catenary degree of 7 , denoted 
0 ( 7 ), is the smallest N such that there exists an A^-chain between any two factorizations of 
7- 

Example 4.2. Consider the numerical semigroup F = (11,36,39). The left-hand picture in 
Figure[^depicts Z(450) and all pairwise distances. Any two factorizations of 450 are connected 
by a 16-chain; one such chain between 2 = ( 6 , 2, 8 ) and w = (24, 3, 2) is depicted with bold 
red edges. Every 16-chain between 2 and w contains the bottom edge, so c(450) = 16. 

The catenary degree can also be computed by examining edges with weight at most 16, 
as depicted in the right-hand picture of Figure [Tj Removing edges labeled 8 and 12 yields a 
(minimal total weight) spanning tree, at which point Lemma 4.5 implies c(450) = 16. 


The graph theoretic description of the catenary degree evident in Example |4.2| is one of 


the core structures used in Algorithm 4.11, We now make this definition explicit. 


with 


Definition 4.3. Given 7 G F, the catenary graph of 7 is the complete graph 
V{K^) = Z(j) whose edges are labeled by the distance function d(—, —). 

To compute the catenary degree of an element 7 G F, one typically first computes Z( 7 ), 
and compiles a list of edges in K-y sorted by distance. From this list of edges, a minimal 
connected graph is constructed one edge at a time, starting with the lowest weight (one can 
use Kruskal’s Algorithm [26] for this). The final edge weight equals the catenary degree. One 


of the primary bottlenecks in this procedure is listing and sorting the edges of since for 
large 7 most edges have weight higher than the catenary degree and thus are never used. 


The main goal of Algorithm 4.11 is to reduce the number of edges by remembering which 
edges were redundant for smaller elements. This is achieved by storing, for each 7 G F, a 
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Figure 1. The catenary graph in F = (11,36,39) from Example 4.2 
Both graphics were created using the computer algebra system SAGE [55] . 


spanning tree of the graph with minimum weight (Definition |4.4[ ), from which the catenary 
degree can be easily recovered (Lemma 4.5). Indeed, the number of edges of is quadratic 
in \ V{K^)\, whereas the number of edges of a spanning tree is linear in \ V{K^)\. 


Definition 4.4. Let G be an undirected graph with integer edge labels. We say a chain 
q = vq,vi, ... ,Vr G y{G) has minimum weight if q has minimum weight among all chains 
in G connecting vq and Vr- A spanning tree T of G has minimum weight if the unique path 
between any two vertices vq and Vr in T is a minimum weight path in G. 


Lemma 4.5. If T is a minimum weight spanning tree of the eatenary graph o/ 7 G F, then 
0(7) equals the highest edge weight in T. 

Proof. Immediate upon unraveling definitions. □ 

The term “minimum weight spanning tree” usually refers to a spanning tree that minimizes 
edge weight sums. Proposition |4.6| verifies that our definition is equivalent to the usual one. 

Proposition 4.6. Fix a graph G, a spanning tree T <Z G with E(T) = {ei,... ,6 ^}, and a 
weight function w on the edges of G. If w{ei) + • • • + w{er) is minimal among all spanning 
trees of G, then T is a minimum weight spanning tree of G. 

Proof. Let p denote a minimum weight chain vo,vi,... ,Vr G I^(G), and let q denote the 
unique path in T from vq to Vr. For i < r, let qi denote the unique path in T connecting 
Vi-i and Vi. Concatenating the chains qi,q 2 , ■ ■ ■ ,Qr produces a chain q' in T connecting vq 
and Vr. Since q is the unique path in T connecting vq and Vr, every edge in q must also be in 
q'. In particular, the edge e of highest weight in q must lie in one of the chains qm. As such, 
if w{p) < w{q), then w{vm-i,Vm) < w{p) < w{q) = w{e), so deleting e from T and adding 
{vm-i,Vm) produces a spanning tree of lesser weight. This completes the proof. □ 

Now that we have identified minimum weight spanning trees as sufficient for computing 
the catenary degree, it remains to see how a minimum weight spanning tree for can be 
constructed from minimum weight spanning trees of its divisors. We accomplish this task 
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using cover morphisms (Definition |4.7| ), which played a crucial role in the dynamic algorithms 
presented in [3] . The key observation is that cover morphisms also preserve edge weights when 
applied to adjacent factorizations in (Remark |4.8[). 


Definition 4.7. A cover morphism is a map ipi : K-^ that acts on the vertex set 

V{K^-ai) = Z (7 — ai) by incrementing the i-th coordinate of each factorization. 


Remark 4.8. When applied to a pair (z, w) of factorizations in the cover morphism 

'ij^i applies the same translation to both z and w, so the resulting factorizations {z', w') satisfy 

z' — gcd( 2 ;^, = z — gcd( 2 ;, vS) and w — gcd( 2 ;^, vJ^ = w — gcd(z, w). 

In particular, d{z\ w') = d(z, w), so preserves edge weights (stated another way, ijji defines 
an isometric embedding of K-y^ai into K^). Additionally, an edge {z',w') G E{K^) lies in the 
image of precisely when i G supp( 2 :') Pi supp(t(;'). Here, supp(z) denotes the support of z, 
that is, the set of indices of the nonzero coordinates of z. 


As a final step. Theorem 4. 10| implies that in addition to edges arising from the minimum 
weight spanning trees of K^^ai, ■ ■ ■, ^ 7 - 0 *, j some supplementary edges from the Graver basis 
of r are needed to obtain a minimum weight spanning tree for K^. The edges from the Graver 
basis of r can be thought of as a base case for the dynamic process, in the sense that every edge 
used in a minimum weight spanning tree is a translation of a Graver basis element. Although 
the Graver basis of T is computationally expensive (albeit necessary; see Example 4.9), it 
only needs to be calculated once for the monoid. 


Example 4.9. Although a minimal presentation for T produces enough edges via translation 
to connect K-y for every 7 G E, these edges are not sufficient for computing catenary degrees. 
Eor example, in the numerical semigroup E = (11,36,39) from Example 4.2, the only edge 
e = ((9, 7,0), (0,0, 9)) in 7^351 with weight c(351) = 16 does not appear in any minimal 
presentation. In fact, the weight 16 edge in Eigure[^is a translation of e. It is this phenomenon 
that forces Algorithm |4.11| to rely on a Graver basis for E. 


Theorem 4.10. For each i <k, let Ti C Ky denote the image in Ky of a minimum weight 
spanning tree of Ky-cti under the cover morphism : Ky-^- —)• Ky. Additionally, let 

E' = R({(z, w) : zioi H-h ZfcOfc = wiai ^ -h Wkak}) n E{Ky) 

he the elements of the Graver basis o/E that lie in Z{'y) x Z(7). The graph 

G' = (Z(7),E(ri)U---UE(rfc)UE') 

is connected, and any minimum weight spanning tree of G' is minimum weight in Ky. 


Proof. Eix z,w ^ ^(7). If j G supp(z) n supp(u;), then a path in Tj connects z and w. 
Alternatively, if z and w have disjoint support and {z,w) ^ E' , then {z,w) can be expressed 
as a sum (z, re) = (zi, z() + • • • + (z^, z(.) of Graver basis elements, so the chain q given by 

Zl + Z2 + • • • + Zr-, Z]^ + Z2 + • • • + Zj., ■ ■ • ; Z]^ + Z 2 + • • • + Z,. 

connects z and w in Ky. Since each edge in q connects factorizations with common support, 
z and w are connected in G'. 

Next, let T’ be any minimum weight spanning tree of G' . It suffices to prove that 
the unique path q' in T' connecting {z,w) G E{Ky) \ E{T') has weight at most d(z,u;). 
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If j G supp(^) nsupp(t(;), then the unique path from z to te in Tj has weight at most d( 2 ;, w), 
and thus so does q'. On the other hand, if z and w have disjoint support, then the i-th edge 
in the chain q constructed in the above paragraph has weight 

d(zj, z') = max(|zi|, \z[\) < max(|z|, \w\) = d(z, w), 
meaning q' has weight at most max{d(zj, z') : i < r} < d(z, w). This completes the proof. □ 


We are now ready to state Algorithm 4.11[ 


Algorithm 4.11. Finds 0 ( 7 ) for any affine semigroup element 7 G F = (oi,..., Uk) C 
function CatenaryDegreeOfAffineSemigroupElement(7, F) 

T ^ MinimumWeightSpanningTree(7, F) 
return max{d(z, ?/;) : {z,w) G E{T)}. 

end function 


H ^ H({(z, w) : zioi H-h ZkUk = wiai H-h tCfcOfc}) 

function MinimumWeightSpanningTree(7, F) 
if Ty not already computed then 
if |Z( 7 )| < 1 then 

r^^(z(7),0) 

else 

E' ^ if n E{K^) 

for all j = 1,..., A: do 

Ti ^ •0i(MlNIMUMWEIGHTSPANNINGTREE(7 - ai,F)) 

end for 

T-y ^ Kruskal(Z( 7), F;(ri) U • • • U E { Tk ) U E ') 

end if 
end if 
return 
end function 


Theorem 4.12. For any ajfine semigroup element 7 G F, Algorithm 4-11 returns 0 ( 7 ). 


Proof. Proposition 4.6 and Theorem 4.10 ensure that MinimumWeightSpanningTree returns 
a minimum weight spanning tree T C K.y, and by Lemma 4.5 the maximum weight among 
the edges of T equals 0 ( 7 ). □ 


4.1. Implementation and benchmarks. Algorithm 4.11 has now been implemented as 
CatenaryDegreeElementListWRTNumericalSemigroup in the numericalsgps m GAP m 
package. The function definition mirrors that of the dynamic length sets and w-primality 
functions [3], which motivated the results in this section. A version for affine semigroups is 
forthcoming, pending discussion of the function definition. The code is accessible in HU or 
from the official GAP web page http://www.gap-system.org, The development version of 
numericalsgps can be found in https://bitbucket.org/gap-system/numericalsgps, 

We give a few implementation notes before concluding the section with an example. 
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1 . 


2 . 

3. 


If r = (ai,... ,afc) C N is a numerical semigroup, then it is natural to begin with 
7 = 0 and compute 0 ( 7 ) in order. This allows the spanning trees to be stored in a 
ring buffer of size ak (the largest generator), since 7 — is the smallest element that 
needs to be recalled to compute 0 ( 7 ). 

By storing the list of edges in E(T^) ordered by weight and computing the union 
E{Ti)U- ■ ■UE{Tk)UE' via a merge of sorted lists, the first (and most computationally- 
intensive) step in Kruskal’s Algorithm runs in linear time. 

At no point must Z('j) be explicitly computed. The base case 7 ( 7 ) < 1 consists of 
precisely those elements for which E{Ti), ..., E{Tk), and E' are all empty, and in 
all other cases, 7 ( 7 ) is simply the set of edge endpoints by Theorem 4.101 This is 


essentially [H Algorithm 3.3] for dynamically computing sets of factorizations. 


Example 4.13. Given here are runtimes for three separate methods of computing 0 ( 7 ) in 
r = (11, 23, 27, 31,43) C N for 7 < 500. 


gap> s:=NumericalSemigroup(ll,23,27,31,43);; 
gap> c:=CatenaryDegreeOfNumericalSemigroup(s);; 
gap> 1:=Intersection([1..500], s);; 

gap> lc:=List(l, x->CatenaryDegreeOfElementInNumericalSemigroup(x,s));; time; 
270390 

gap> ln:=List(l, x->CatenaryDegreeOfElement(x,s,c));; time; 

29380 

gap> ld:=CatenaryDegreeElementListWRTNumericalSemigroup(l,s);; time; 

6630 


The function CatenaryDegreeOf ElementInNumericalSemigroup (original implementation 
in the package) computes 0 ( 7 ) by computing the complete graph K-y and then removing edges 
by weight, whereas the function CatenaryDegreeOfElement constructs a spanning tree for 
Ky using Kruskal’s Algorithm (as discussed immediately following Definition 4.3). The lat¬ 
ter will replace the implementation of the former in the next release of the numericalsgps 
package. Although CatenaryDegreeOfElement is an improvement over the original imple¬ 
mentation, the dynamic algorithm is still considerably faster for large 7 . 


5 . The tame degree 


A variant of the catenary degree, the tame degree also uses distance (as in Definition |4.1[ ) to 
measure how wild the factorizations of monoid elements are [22]. In this section, we present 
Algorithm |5.10 an improved method for computing the tame degree of a full affine semigroup 
(Definition 5.5). We begin by presenting Theorem 5.2, the main new theoretical result used 
in Algorithm 5.10, followed by a discussion of the algorithm. One of the primary motivations 
for having a dedicated algorithm for this family of affine semigroups is to examine the tame 
degree of block monoids; in Examples [5.15[|5 .17 we apply our algorithm in this setting. 


Definition 5.1. Fix T = (oi,..., Ofc) C N^. Given 7 G T and i < k such that 7 — Oj G T, 
the tame degree of 7 with respect to Oj, denoted tj(7), is defined as the least G N U {00} 
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such that for every z G Z( 7 ), there exists w G 7 ( 7 ) with Wi > 0 and 6{z,w) < N. The tame 
degree of T with respect to ai is then defined as 

ti(r) = sup{ti( 7 ) : 7 - a* G r}. 

Lastly, the tame degree of T is dehned as t(r) = max{ti(r): i < k}. 

In what follows, we denote by ‘<’ the usual partial ordering in that is, x < y if 
y — X G N^, or, equivalently, if Xi < yi for all i < k. Notice that x < y is also equivalent to 
yGx + N*^ = {x + 2 :: z € N^}. 

Theorem 5.2. Let T = (oi,..., ak) C and fix i < k. We have 

ti(r) = sup{ti( 7 ) : 7 G Mi}, 
where Mi = {(/9r(a^) : x G Minimals<Z(ai + T)}. 


Proof. Since Mj C T, clearly sup{tj(m) : m G Mj} < t(r). Conversely, fix 7 G a* + T 
and 2 ; G Z(7). If Zi > 0, then choose w = z in the definition of tame degree. Otherwise, 
2 ; G Z(aj + r) since 7 G a* + T, and consequently, there exists z' G Minimals<Z(ai + T) such 
that z' < z. Set 7' = (pr{z'). By the definition of tj(— ), there exists w' G Z{y) n (ai + N^) 
with d(2;', w') < ^(7'). Let w = w' + (z — z'). Upon verifying that 

+ Pr(z - z') = ipr{z') + y:=r{z - z') = (pr{z) = 7 , 

d( 2 :,ui) = d( 2 ;',ui') < tj( 7 '), and w G Z{'y) n (a* + N*^), it follows that tj(n) < tj(m). This 
implies tj( 7 ) < sup{tj( 7 ') : 7 ' G Mi}, and consequently tj(r) < sup{tj( 7 ') : 7 ' G Mi}. □ 


Remark 5.3. Theorem |5.2| holds in more generality with a nearly identical proof; the only 
requirement is that T be atomic. We specialize to the affine setting here to simplify notation 
since Algorithm 5.10 requires this assumption; see [22] for the general definitions. 


Remark 5.4. Comparing Theorem 5.2 to [H Proposition 3.3], it is not a coincidence that the 
elements needed to compute the tame degree of T are the same as those needed to compute the 
w-primality invariant (see [U Dehnition 3.1]). Some evidence of this connection was already 
observed for numerical semigroups; see, for instance, I?] Corollary 3] and [U Remark 5.9]. 


Arranging the atoms ai,... ,ak of T as the columns of a matrix A, the set Z( 7 ) of factor¬ 
izations of 7 G T coincides with the nonnegative integer solutions to the linear system 

Ax = 7 . 

Prior to Theorem |5.2[ the tame degree of P was typically computed by first computing the 
Graver basis of A |9|. We propose here an alternative method, instead computing the sets 

Mi = Minimals<Z(ai -|- P) and Mi = {(py(z) : z G Mi} 

for each i G {1,..., k}. By Dickson’s lemma, both Mi and Mi have finitely many elements. 

For the calculation of Mi one can use [321 Algorithm 16] (this is precisely how w-primality 
is computed in m). However, it turns out that Normaliz mu performs this task faster. 
To obtain Mi, we hrst compute the minimal nonnegative integer solutions of 

(-4 I - A ) (}}} = o., 
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and then project on the first k coordinates (here, {A \ —A) is the matrix having as columns 
the columns of A followed by those of —A). From the resulting (finite) set, we simply take 
those elements that are minimal with respect to <. 

Since the Graver basis of A coincides with the minimal nonnegative integer solutions of 

0 ^ 

it would seem that our approach has no signihcant computational advantage over the proce¬ 
dure given in [9]. However, combining with [H Corollary 3.5] (stated here as Proposition |5.7| ) 
yields signihcant performance improvements for full affine semigroups. 

Definition 5.5. An affine semigroup F C N'^ is full if G(F) n = F, where G(F) denotes 
the subgroup of generated by F. 


Remark 5.6. Up to isomorphism, every full affine semigroup is a (reduced hnitely generated) 
Krull monoid [251 Proposition 2]. Also, the class of reduced Krull monoids coincides with the 
class of saturated submonoids of free monoids [221 Theorem 2.4.8]. 

Proposition 5.7 (jH Corollary 3.5]). //F = (oi,..., Ofc) C is full, then for every 7 G F, 

Z (7 -|- F) = {x G : 7 < iprix)}- 


Example 5.8. The “full” hypothesis in Proposition |5.7| cannot be omitted. For example, 
the numerical semigroup F = (3, 5, 7) C N is not full since 5 — 3 = 2 G (G(F) n N) \ F (in fact, 
the only full numerical semigroup is N). In this example, we see that 

Minimals<Z(3 + F) = {(0, 0, 2), (0,1,1), (0, 2, 0), (1,0,0)}, 
even though Minimals<{x G : 3xi -|- 5 x 2 + 7 x 3 > 3} = {(1, 0,0), (0,1, 0), (0, 0,1)}. 


Recalling A = (oi | • • • | Ofc) from above, Proposition 5.7 states that 

Z(7 + F) = {x G : Ax > 7}, 

meaning that can be computed as a Hilbert basis. At this point, in order to compute tj(F) 
using Theorem 5.2 it remains to compute tj(7) for each 7 G Mj. To this end. Proposition ]5.9| 
implies that for hxed z G Minimals<Z(aj-|-F) with Zi = 0 , we can hnd the closest factorization 
rc G to z in the hber of Pt{z) with Wi 0 hy minimizing |r(;| subject to the constraints 
Aw = Az and w ■ z = 0. In particular, this expresses tj( 7 ) as the result of an integer linear 
programming problem. 


Proposition 5.9. Let F = {a\, ..., a^) C fix i < k, and fix z £ Minimals<Z(ai -|- F). 
If Zi = 0 and w £ fulfills Wi 0 and Aw = Az, then z-w = 0 and d(z,t(;) = max{|z|, |rc|}. 


Proof Assume to the contrary that u = gcd(z, w) 7 ^ 0. As Zi = 0, we deduce that Ui = 0. 
This implies that the i-th coordinate of w — u is nonzero, and consequently w — u £ Z(aj-|-F). 
Hence z — u £ Z(ai + F) and z — u < z, contradicting the minimality of z. □ 


The last necessary observation is that e* G Minimals<Z(aj -|-F), where e* is the Ah column 
vector of the k x k identity matrix, so it suffices to consider those 2; G Minimals<Z(Q;j -|- F) 


satisfying z* = 0 when computing tj(F). We now summarize our procedure in Algorithm 5.10 
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Algorithm 5.10. Computes ti(r) for a full affine semigroup F = (ai,..., ak) C N'^. 
function TAMEDEGREEOFFULLAFFINESEMIGROUP(r,i) 

H ^ minimal solutions z of Ax > ai with respect to < satisfying Zi = 0 
if H is empty then 

return 0 (and thus stop) 
end if 
P ^ 0 

for all z ^ H do 

minimize yi + ■ ■ ■ + Vk subject to Ay = Az and y • z = 0 
P^PVJ{y} 

end for 

return max{|z|, |y| : z G iF,y G P}. 

end function 


Theorem 5.11. For any full affine semigroup F, Algorithm 5.10 returns tj(F). 


Proof. Theorem 5.2, Propositions 5.7 and 5.9, and the discussion in between. 


□ 


Remark 5.12. Notice that if all solutions y obtained in Algorithm 5.10 satisfy |y| < \z\, 
then tj(F) will coincide with the w-primality of Oj, u}{ai) (see [U Proposition 3.3]). Indeed, 
this provides an alternative proof that io{ai) < tj(F). 


5.1. Implementation and benchmarks. We have implemented Algorithm 5.10 in the 
numericalsgps [14j GAP [T7] package. The function TameDegreeOfAffineSemigroup therein 
has two main methods: one for affine semigroups in general, and another for full affine semi¬ 
groups (affine semigroups with the attribute HasEquations). The code is accessible in [T3] 
or from the official GAP web page http://www.gap-systeni.org, The development version 
of numericalsgps can be found in https://bitbucket.org/gap-system/numericalsgps, 
We conclude this section by applying our implementation of Algorithm 5.10 in several 
examples of block monoids (Definition |5.13[) . 


Definition 5.13. Fix an abelian group (G, -t-) and a finite subset H = {yi,... ,y„} C G\{0}. 
The block monoid associated to H is the affine semigroup 

B{H) = {(xi,..., x„) G : xiyi + • • • + XnQn = 0} C N^. 

Denote by B{G) the block monoid of the set of nonzero elements of G. 


Remark 5.14. It is easy to show that B{H) is a full affine semigroup. Many factorization 
properties of a monoid can be derived from the factorization properties of the block monoid 
associated to its class group (see [22] for a thorough treatment). 

Example 5.15. We begin by computing the tame degree of F = BifL^. If B denotes the 
matrix with columns the nonzero elements of Z^, then F is the submonoid of of nonnegative 
integer solutions of the system 

Bx = 0 mod 2. 

The tame degree of F can be computed in GAP as follows. 
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gap> m: = [ [ 0 , 0 , 1 ] , [ 0 , 1 , 0 ], [ 0 , 1 , 1 ], [ 1 , 0 , 0 ], 

[ 1 , 0 , 1 ], [ 1 , 1 , 0 ], [ 1 , 1 , 1 ]]; 
gap> a:=AffineSemigroupC'equations",[TransposedMat(m),[2,2,2]]); 
<Affine semigroup> 


gap> at:=GeneratorsOfAffineSemigroupCa); 


[ [ 

0 , 

0 , 

0 , 

0 , 

0 , 

2 , 

0 

]. 

[ 

0 , 

0 , 

0 , 

0 , 

0 , 

0 , 

2 

], 

[ 

0 , 

0 , 

0 , 

0 , 

2 , 

0 , 

0 

], 

[ 

0 , 

0 , 

0 , 

2 , 

0 , 

0 , 

0 

]. 

[ 

0 , 

0 , 

2 , 

0 , 

0 , 

0 , 

0 

], 

[ 

0 , 

2 , 

0 , 

0 , 

0 , 

0 , 

0 

], 

[ 

2 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 

]. 

[ 

0 , 

0 , 

1 , 

0 , 

1 , 

1 , 

0 

], 

[ 

0 , 

0 , 

1 , 

1 , 

0 , 

0 , 

1 

], 

[ 

0 , 

1 , 

0 , 

0 , 

1 , 

0 , 

1 

]. 

[ 

0 , 

1 , 

0 , 

1 , 

0 , 

1 , 

0 

], 

[ 

1 , 

0 , 

0 , 

0 , 

0 , 

1 , 

1 

], 

[ 

1 , 

0 , 

0 , 

1 , 

1 , 

0 , 

0 

]. 

[ 

1 , 

1 , 

1 , 

0 , 

0 , 

0 , 

0 

], 

[ 

0 , 

0 , 

0 , 

1 , 

1 , 

1 , 

1 

], 

[ 

0 , 

1 , 

1 , 

0 , 

0 , 

1 , 

1 

]. 

[ 

0 , 

1 , 

1 , 

1 , 

1 , 

0 , 

0 

], 

[ 

1 , 

0 , 

1 , 

0 , 

1 , 

0 , 

1 

], 

[ 

1 , 

0 , 

1 , 

1 , 

0 , 

1 , 

0 

]. 

[ 

1 , 

1 , 

0 , 

0 , 

1 , 

1 , 

0 

], 

[ 

1 , 

1 , 

0 , 

1 , 

0 , 

0 , 

1 

] ] 


gap> TameDegreeOfAffineSemigroup(a); 
4 


We use the package Normalizlnterf ace |23j that utilizes the Normaliz library to compute 
the Hilbert bases. Although 4ti2 [1] can also be used to compute Graver and Hilbert bases via 
the 4ti2gap [2] package, in this example 4ti2 was slower than Normaliz. The computation 
of t(r) took 4955 ms on an i7 laptop with 16 GB of memory (Normaliz was compiled without 
OpenMP, and thus was running in a single thread). For comparison, an implementation of the 
procedure described in [9] took 1,048,757 ms. 


The runtime can actually be improved further. 
t(r) it suffices to compute ti(r) with \ai\ = 4. 
minimalElementsPrincipalldealOfFullAff ineSemigroup implements Proposition 5.7 


According to m, in order to compute 
This is accomplished as follows, where 


gap> u:=First(at,x->Sum(x)=4); 

[ 0 , 0 , 0 , 1 , 1 , 1 , 1 ] 

gap> Mu:=minimalElementsPrincipalIdealOfFullAffineSemigroup(u,a);; 
gap> facts:=List(Mu, x->FactorizationsVectorWRTList(x,at));; 
gap> Set(facts,TameDegreeOfSetOfFactorizations); 

[ 0, 2, 3, 4 ] 

It took 2 ms to compute at (the set of atoms of T), and 12 ms to compute Mu. It took 
another 420 ms to complete the third line above, and the final line finished in 425 ms. In 
total, it takes less than a second to compute the tame degree of T. 


Example 5.16. Running Algorithm 5.10 for T = H(Z 2 x Z 3 ) yields t(r) = 8 after 19214 ms; 
execution of the procedure presented in [9] was killed after several hours. 


Example 5.17. As the size of G increases, the number of atoms (i.e. the ambient dimension 
of B{G)) quickly makes use of Algorithm 5.10 infeasible. Indeed, the block monoid ^(Zl) has 
323 atoms, and ^(Zl) has 20367 atoms. Although Algorithm 5.10 is considerably faster than 
the general-purpose algorithm introduced in [9], it is still not sufficient to compute t(H(Z|)). 
Given below is the analysis for t(H(Z|)) by using the supercomputer alhambra.ugr.es, 


(1) First attempt: 16 cores with 30 GB of internal memory. 
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Internal memory is quickly exhausted and the system begins using swap memory, decaying 
the number of processors and the user CPU usage. The process is killed before it hnishes. 
(2) Second attempt: 32 cores with 256 GB of internal memory. 


1 " 



After 72 hours of execution, Mu is still not computed. We contacted W. Bruns, B. Ichim 
and C. Soger (developers of Normal iz) and were told (i) the input has too many variables 
and few equations to finish in a reasonable time, and (ii) the cone has tens of millions of 
extremal rays and thus requires extensive memory (> 50 GB) just to store. 

The point of this example is that even though we may consider Z| a “small” group, 
the monoid . 6 (^ 2 ) has 20367 atoms, and this translates to working in an ambient space of 
dimension 323! Thus, Algorithm |5.10 will fail to compute the tame degree for this monoid 
with the computational tools at hand. 
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